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Executive  Summary 


We  describe  the  development  of  PRISM  (Parameterized  Peal-time  /onospheric 
Specification  Model),  a  global  near  real-time  ionospheric  specification  model  intended  for  use  at 
the  Air  Force  Space  Forecast  Center  (AFSFC),  also  known  as  the  50th  Weather  Squadron 
(50  WS).  We  also  describe  PIM  (Parameterized  /onospheric  Model),  the  base  ionospheric  model 
which  forms  the  core  of  PRISM.  PIM  consists  of  a  semi-analytic  representation  of  a 
parameterization  of  four  separate  physically  based  computational  models  of  the  ionosphere. 
PRISM  uses  both  ground  based  and  space  based  data  to  modify  or  update  the  PIM  ion  density 
profiles  to  produce  a  data  driven  specification  of  the  state  of  the  ionosphere.  At  AFSFC  PRISM 
will  provide  hourly  ionospheric  specifications  in  near  real  time.  PRISM  is  capable  of  accepting 
and  using  bottomside  electron  density  profiles  produced  from  automated  true  height  analysis  of 
digital  ionograms,  TEC  measurements  from  GPS  receivers,  satellite  based  in  situ  plasma 
measurements  (electron  density,  ion  composition,  electron  and  ion  temperature,  and  ion  drift 
velocity),  satellite  based  in  situ  observations  of  precipitating  electrons  and  ions  in  the  high  latitude 
regions,  as  well  as  ionospheric  parameters  derived  from  observations  of  airglow  and  auroral 
optical  emissions.  We  describe  the  construction  of  PIM,  and  the  real  time  adjustment  process  of 
PRISM,  as  implemented  in  Version  1.5.  This  is  the  version  that  will  become  operational  at 
AFSFC  in  the  fall  of  1995. 
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1.  INTRODUCTION 


We  have  developed  a  Parameterized  Real-time  Ionospheric  Specification  Model  (PRISM) 
for  the  Air  Force  Air  Weather  Service  (AFAWS)  for  use  at  the  Air  Force  Space  Forecast  Center 
(AFSFC,  also  known  as  the  50th  Weather  Squadron).  The  model  uses  both  ground  based  and 
space  based  data  available  in  near  real  time  to  modify  a  Parameterized  Ionospheric  Model  (PIM) 
and  thus  provide  a  near  real-time  specification  of  the  ionosphere.  PIM  is  a  composite  of  diumally 
reproducible  runs  of  several  physical  ionospheric  models:  (1)  the  Time  Dependent  Ionospheric 
Model  (TDIM)  of  Utah  State  University  (USU)  [Schunk,  1988],  (2)  the  low  latitude  F-region 
model  (LOWLAT)  developed  by  Anderson  [1973],  (3)  the  midlatitude  version  of  LOWLAT 
(called  MIDLAT)  developed  by  D.  N.  Anderson  and  modified  by  D.  T.  Decker,  and  (4)  an  E- 
region  local  chemistry  code  developed  by  D.  T.  Decker  and  incorporating  photoelectrons  using 
the  continuous  slowing  down  method  [. /asperse ,  1982], 

Both  PIM  and  PRISM  can  produce  either  regional  or  global  output.  The  output  grid  of 
latitude  and  longitude  (which  may  be  either  geographic  or  geomagnetic)  is  user  selectable.  At 
AFSFC  the  standard  PRISM  output  is  a  global  grid  with  1°  latitude  spacing  and  5°  longitude 
spacing.  Profile  parameters  (faF2,  hmF2,  TEC,  etc.),  or  complete  electron  density  profiles,  or  both 
may  be  output  to  the  specified  grid. 

PRISM  has  been  delivered  to  the  AFAWS  and  is  presently  undergoing  transition  to 
operational  code.  PIM  is  in  the  public  domain  and  is  available  for  distribution  to  the  ionospheric 
community.  PRISM  is  also  available  to  the  research  community  upon  application  to  the  AFAWS. 
PIM  has  been  distributed  to  about  60  users  worldwide.  PRISM  is  in  use  by  three  research  groups 
in  the  United  States. 
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1.1  Objectives 

Our  primary  objective  was  the  development  of  an  algorithm  for  using  near  real  time 
satellite  and  ground  based  data  to  provide  a  near  real  time  specification  of  the  global  ionosphere. 
The  data  to  be  used  include 

(1)  bottomside  digital  soundings  from  the  AFAWS  Digital  Ionospheric  Sounding  System  (DISS), 

(2)  Total  Electron  Content  (TEC)  data  from  the  AFAWS  Transionospheric  Sensing  System  (TISS), 

(3)  in  situ  plasma  data  (densities,  temperatures,  and  drift  velocities)  from  the  SSIES  instrument  on 
DMSP  satellites, 

(4)  auroral  electron  and  ion  fluxes  from  the  SSI/4  instrument  on  DMSP  satellites,  and 

(5)  electron  density  profile  information  deduced  from  observations  of  airglow  and  auroral  optical 
emissions  by  instruments  (SSUSI  and  SSULI)  expected  to  be  flown  on  future  DMSP  satellites. 

The  need  for  a  global  specification  of  the  state  of  the  ionosphere  is  twofold.  First,  there 
are  operational  systems  that  need  to  correct  for  ionospheric  effects  in  real  time,  or  that  have 
operational  parameters  that  are  affected  by  the  ionosphere  and  must  be  adjusted  in  near  real  time. 
Second,  the  operation  of  many  systems  could  be  optimized  if  accurate  forecasts  of  ionospheric 
conditions  were  available  because  this  allows  the  operational  parameters  to  be  chosen  ahead  of 
time.  Any  ionospheric  forecast  algorithm  will  require  an  accurate  specification  of  the  current  state 
of  the  ionosphere  as  an  initial  condition,  which  PRISM  provides. 

In  addition  to  real  time  needs,  many  system  operators  need  post-event  analysis  to 
determine  whether  operational  problems  or  outages  were  caused  by  system  problems  or  by 
environmental  conditions.  PRISM  will  be  used  for  this  purpose  at  AFSFC  as  well. 


1.2  Approach 

Ideally,  the  specification  of  the  current  state  of  the  ionosphere  would  be  obtained  directly 
from  real  time  observations  from  a  dense  network  of  satellite  and  ground  based  instruments. 
Unfortunately,  the  complexity  and  spatial  extent  of  the  ionosphere  precludes  the  deployment  of  a 
sufficiently  dense  network  of  observing  instruments.  Therefore,  any  practical  ionospheric 
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specification  algorithm  must  be  based  on  an  ionospheric  model  with  parameters  that  can  be 
adjusted  on  the  basis  of  near  real  time  data.  Two  approaches  are  possible:  (1)  statistical  or 
empirical  climatological  models  or  (2)  numerical  simulations  based  on  physical  models.  (In  this 
paper  “physical  model”  is  synonymous  with  “numerical  model.”)  For  reasons  described  below, 
we  have  chosen  the  second  approach  (physical  models).  However,  practical  considerations 
(primarily  computational  speed)  dictate  that  the  algorithms  implemented  at  the  Space  Forecast 
Center  be  based  on  parameterized  versions  of  the  physical  models  (“theoretical  climatology”). 

We  feel  strongly  that  a  comprehensive  physical  model  of  ionospheric  processes  can 
produce  more  accurate  specifications  and  forecasts  than  can  statistical  or  climatological  models. 
The  causal  relationship  between  easily  monitored  solar  and  geophysical  parameters  (e.g.,  Kp  F101, 
etc.)  and  a  particular  ionospheric  configuration  is  very  complex.  Any  organization  of  historical 
ionospheric  data  inevitably  averages  over  a  variety  of  ionospheric  configurations  corresponding  to 
similar  values  of  the  chosen  set  of  solar-geophysical  parameters  (usually  only  one  or  two).  The 
result  is  that  spatial  structure  tends  to  be  smeared  out  or  smoothed  over,  and  the  resulting  model 
is  unrepresentative  of  the  instantaneous  ionosphere.  If  a  physical  model  contains  all  of  the 
relevant  physics,  and  if  the  inputs  are  realistic,  then  it  will  produce  more  realistic  representations 
of  instantaneous  ionospheric  structure.  However,  there  is  a  difference  between  a  realistic 
representation  and  an  accurate  one. 

In  order  to  accurately  simulate  a  time  dependent  phenomenon  like  the  ionosphere,  a 
physical  model  needs  an  accurate  specification  of  the  initial  conditions  and  an  accurate 
representation  of  the  energy  and  momentum  flux  at  the  boundaries.  For  the  purposes  of  providing 
a  specification  model,  it  is  the  energy  and  momentum  input  that  is  crucial.  If  the  model  is  run  long 
enough,  the  effects  of  the  initial  conditions  are  lost  and  the  present  state  of  the  model  depends 
only  on  the  recent  history  of  the  energy  and  momentum  input.  These  include  the  solar  EUV  (the 
primary  source  of  ionization  outside  the  auroral  zone),  high  latitude  heating  of  the  thermosphere 
(which  affects  the  global  circulation  of  the  thermosphere),  high  latitude  convection,  and  low 
latitude  dynamo  electric  fields.  While  the  temporal  and  spatial  resolution  of  the  observations  of 
these  quantities  are  expected  to  improve  in  the  future,  they  will  probably  never  be  sufficient  to 
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allow  accurate  simulation  of  the  ionosphere  without  additional  data.  As  a  practical  matter, 
ionospheric  simulations  must  be,  and  will  remain,  iterative  in  nature:  The  energy  and  momentum 
input  parameters  are  adjusted  until  the  simulation  agrees  with  observations  of  ionospheric 
parameters  to  some  level  of  accuracy. 

A  practical  consequence  of  this  situation  is  that  the  production  of  an  accurate  ionospheric 
specification  based  on  a  numerical  simulation  requires  that  the  physical  model  be  run  several 
(perhaps  many)  times  until  its  output  agrees  with  the  available  observations.  The  requirement  that 
the  model  cover  a  sufficient  time  period  to  allow  the  transient  effects  of  the  initial  conditions  to 
damp  out  implies  that  the  physical  model  must  run  much  faster  than  the  system  it  is  simulating. 
At  the  present  time,  given  practical  (i.e.,  financial)  limits  on  the  available  computing  power,  this  is 
not  possible.  Consequently,  we  have  adopted  a  modified  approach  in  which  the  physical  models 
are  parameterized  in  terms  of  solar  and  geophysical  parameters.  It  is  these  parameterized  models 
rather  than  the  original  physical  models  which  are  to  be  adjusted  according  to  the  real  time 
ionospheric  data. 

There  is  a  superficial  similarity  between  our  approach  and  a  climatological  approach.  The 
difference,  however,  is  that  we  begin  with  a  more  realistic  representation  of  the  spatial  structure 
of  the  ionosphere  than  climatological  models  can  provide.  The  parameter  adjustment  process 
should  not  compromise  this  advantage.  In  the  future,  as  more  powerful  computers  and  more 
efficient  model  algorithms  become  available,  the  parameterized  models  can  be  replaced  by  actual 
physical  models  to  produce  more  accurate  specifications. 
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2.  THE  PHYSICAL  MODELS 


Four  separate  physical  models  were  used  as  the  basis  of  PRISM:  (1)  a  low  latitude  F  layer 
model  (LOWLAT),  (2)  a  midlatitude  F  layer  model  (MIDLAT),  (3)  a  combined  low  and  middle 
latitude  E  layer  model  (ECSD),  and  (4)  a  high  latitude  E  and  F  layer  model  (TDIM).  All  four 
models  are  based  on  a  tilted  dipole  representation  of  the  geomagnetic  field  and  a  corresponding 
magnetic  coordinate  system.  (Hereafter,  “latitude”  means  “magnetic  latitude”  unless  otherwise 
noted.)  All  four  models  use  the  MSIS-86  neutral  atmosphere  model  [Hedin,  1987].  Chemical 
reaction  rates,  collision  frequencies,  and  similar  data  are  consistent  among  all  the  models. 


2.1  The  Low  Latitude  F  Layer  Model 

The  low  latitude  F  region  model  (LOWLAT)  was  originally  developed  by  Anderson, 
[1973].  (See  also  Moffett,  [1979]).  It  solves  the  diffusion  equation  for  0+  along  a  magnetic  flux 
tube.  Normally,  the  entire  flux  tube  is  calculated  with  chemical  equilibrium  boundary  conditions 
at  both  feet  of  the  flux  tube.  A  large  number  of  flux  tubes  must  be  calculated  in  order  to  build  up 
an  altitude  profile. 

Since  heat  transport  is  not  included  in  this  model,  ion  and  electron  temperature  models 
must  be  used.  For  the  PRISM  development  effort  we  chose  the  temperature  model  of  Brace  and 
Theis  [1981].  The  Horizontal  Wind  Model  (HWM)  of  Hedin  [1988]  was  used  to  describe 
thermospheric  winds. 

The  critical  feature  incorporated  in  the  low  latitude  model  is  the  dynamo  electric  field. 
The  horizontal  component  of  this  field  drives  upward  convection  through  ExB  drift,  and  this  can 
significantly  modify  profile  shapes  and  densities.  This  phenomenon  is  responsible  for  the 
equatorial  anomaly,  crests  in  ionization  on  either  side  of  the  magnetic  equator  at  ±15—20° 
magnetic  latitude.  In  the  current  version  of  PRISM  (Version  1.5)  the  ExB  driven  vertical  drift 
used  for  these  calculations  was  based  on  the  empirical  models  derived  from  data  from  the 
Atmospheric  Explorer-E  (AE-E)  satellite  [Fejer  et  al,  1995],  which  are  consistent  with  the  drifts 
measured  at  Jicamarca  [Fejer,  1981;  Fejer  et  al.,  1989]  but  include  longitudinal  variations  as  well. 
We  used  the  Fejer  et  al.  [1995]  empirical  drifts  for  moderate  and  high  solar  activity.  Following 
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the  their  discussion,  we  modified  these  drifts  by  reducing  or  eliminating  the  pre-reversal 
enhancement  for  low  solar  activity.  Horizontal  drifts  were  neglected  in  the  PRISM  runs. 

Since  its  original  development  this  model  has  undergone  extensive  validation  by 
comparison  with  data.  The  most  recent  such  comparison  is  Preble  et  al.  [1994]  using  electron 
density  profiles  measured  by  the  incoherent  scatter  radar  facility  at  Jicamarca,  Peru. 


2.2  The  Midlatitude  F  Layer  Model 

The  midlatitude  F  region  model  (MIDLAT)  is  the  same  as  the  low  latitude  version,  except 
that  the  dynamo  electric  field  is  not  included.  Complete  flux  tubes  are  followed,  but  neither 
horizontal  nor  vertical  convection  is  included.  The  computer  resource  requirements  of  MIDLAT 
are  far  less  than  those  of  LOWLAT.  As  long  as  the  boundary  between  low  and  middle  latitudes  is 
chosen  so  that  the  electric  field  is  negligible  on  the  boundary  flux  tubes,  the  two  models  give 
identical  results  at  the  boundary  ensuring  continuity  across  that  boundary.  For  the  PRISM 
development  effort  we  used  the  same  temperature  model  [Brace  and  Theis,  1981]  and  the  same 
thermospheric  wind  model  [Hedin,  1988].  For  appropriate  production,  loss,  and  diffusion  rates 
for  both  LOWLAT  and  MIDLAT,  see  Decker  et  al.  [1994], 


2.3  The  Low  and  Midlatitude  E  Layer  Model 

The  low  and  mid-  latitude  E  region  model  (ECSD)  was  developed  by  Dwight  T.  Decker 
and  John  R.  Jasperse  and  incorporates  photoelectrons  calculated  using  the  continuous  slowing 
down  (CSD)  approximation  [Jasperse,  1982].  Ion  concentrations  are  calculated  assuming  local 
chemical  equilibrium  A  small  nighttime  source  is  included  to  ensure  that  an  E  layer  is  maintained 
throughout  the  night. 
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2.4  The  High  Latitude  Model 


The  high  latitude  model  (incorporating  both  E  and  F  layers)  is  the  Utah  State  University 
(USU)  Time  Dependent  Ionospheric  Model  (TDIM).  (See  Schunk  [1988]  for  a  review.)  This 
model  is  similar  to  the  low  and  middle  latitude  models  except  that  the  flux  tubes  are  truncated  and 
a  flux  boundary  condition  is  applied  at  the  top.  In  addition,  the  flux  tubes  move  under  the 
influence  of  the  high  latitude  convection  electric  field.  In  the  low  latitudes,  because  the  magnetic 
field  is  mainly  horizontal,  the  effect  of  the  electric  field  is  primarily  to  move  the  ionization  in 
altitude.  In  contrast,  the  high  latitude  magnetic  field  is  mainly  vertical,  and  the  electric  field  driven 
convection  is  horizontal.  Like  LOWLAT,  this  model  has  a  long  history  and  has  been  validated  by 
numerous  comparisons  with  data. 

TDIM  includes  an  E-layer  model  that  incorporates  the  effects  of  ionization  by 
precipitating  auroral  particles.  The  ion  production  rates  used  were  calculated  using  the  B3C 
electron  transport  code  [Strickland,  1976;  Strickland  et  al.,  1994]  and  incident  electron  spectra 
representative  of  DMSP  SSJ/5  data.  The  characteristics  of  the  electron  spectra  were  taken  from 
the  Hardy  et  al.  [1987]  electron  precipitation  model.  The  high  latitude  convection  patterns  were 
those  developed  by  Heppner  and  Maynard  [1987]  for  southward  directed  Bz. 
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3.  PARAMETERIZATION  OF  THE  PHYSICAL  MODELS 


Parameterization  of  the  physical  models  proceeded  in  two  steps.  First,  the  models  were 
used  to  generate  a  number  of  "databases"  for  a  discrete  set  of  geophysical  conditions.  Each 
database  consists  of  ion  density  profiles  on  a  discrete  grid  of  latitudes  and  longitudes  for  a  24  hour 
period  in  UT.  Second,  to  reduce  storage  requirements,  the  databases  were  approximated  with 
semi-analytic  functions.  These  two  processes  are  described  in  the  following  subsections. 

3.1  Geophysical  Parameters 

All  the  physical  models  were  parameterized  in  terms  of  season  and  solar  activity.  The 
middle  and  high  latitude  models  were  also  parameterized  in  terms  of  magnetic  activity,  while  the 
high  latitude  model  was  additionally  parameterized  in  terms  of  the  sign  of  the  interplanetary 
magnetic  field  component  By.  (The  high  latitude  model  was  only  run  using  Bz  southward. 

Northward  Bz  conditions  are  modeled  using  the  low  magnetic  activity  databases.)  For  the  middle 
and  low  latitudes,  the  F  layer  (CL)  and  the  E  layer  (NO+and  Oj)  were  computed  and 
parameterized  separately.  The  high  latitude  model  (TDIM)  produced  all  three  ions 
simultaneously. 

Due  to  time  and  computer  resource  limitations,  only  a  few  values  of  each  parameter  were 
used.  The  season  "values"  are  the  June  and  December  solstices  and  the  March  equinox  (which 
also  "stands  in"  for  the  September  equinox).  We  expect  to  change  from  seasonal  to  monthly 
values  in  the  next  versions  of  PIM  and  PRISM.  The  values  of  the  other  parameters  are 
summarized  for  each  latitude  region  in  Table  1.  Note  that  the  USU  TDIM  and  LOWLAT 
produce  output  in  magnetic  local  time  (MLT),  while  MIDLAT  and  ECSD  produce  output  in 
magnetic  longitude.  Since  the  two  coordinates  are  readily  interconvertible,  we  will  ignore  the 
distinction  and  refer  only  to  magnetic  longitude  in  the  description  that  follows. 
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Tab 


e  1:  Geophysical  Parameter  Values 


Solar  Activity 
(^0.7) 

Magnetic 

Activity 

(Kp) 

IMF  By 

direction 

Number  of 

databases 

Low  Latitude  F  layer 

70, 130, 210 

1 

N/A 

36° 

Midlatitude  Flayer 

70, 130, 210 

■ 

N/A 

546 

Low  &  Midlatitude  E  layer 

70,  130, 210 

N/A 

54c 

High  Latitude  E  &  Flayer 

70, 130, 210 

1,  3.5,  6 

+,  “ 

324 d 

a3  seasons  x  3  solar  activities  x  4 
b3  seasons  x  3  solar  activities  x  3 
c3  seasons  x  3  solar  activities  x  3 
d3  seasons  x  3  solar  activities  x  3 


longitude  sectors 

magnetic  activities  x  2  hemispheres 
magnetic  activities  x  2  species 

magnetic  activitiesx  2  By’sx  3  species  x  2  hemispheres 


3.2  Representation  of  the  Databases 

When  the  models  are  run  for  any  one  set  of  geophysical  parameters  (e.g.,  June,  FW1  -  130, 
Kp  =  1),  they  produce  ion  densities  (0+,  NO+,  and  0+_)  on  a  four  dimensional  grid.  MIDLAT  and 
ECSD  use  a  grid  of  magnetic  latitude  ( X ),  magnetic  longitude  (<p),  altitude(z),  and  Universal 
Time  (x ).  TDIM  uses  magnetic  local  time  (MLT  or  i|/ )  instead  of  magnetic  longitude,  while 
LOWLAT  uses  MLT  instead  of  UT.  In  order  to  make  this  mass  of  numbers  more  manageable, 
we  produced  a  semi-analytical  representation  of  each  database.  The  space  and  time  grid 
parameters  are  summarized  for  each  latitude  region  in  Table2. 

Due  to  the  computer  resource  requirements  of  the  low  latitude  F  layer  code,  it  was  used  to 
generate  databases  at  four  discrete  longitudes  (corresponding  to  longitude  sectors  for  which  ExB 
drift  measurements  were  available).  Each  longitude  sector  was  parameterized  separately,  and  the 
necessary  longitude  interpolation  is  carried  out  in  PIM  and  PRISM  during  execution,  as  described 
below.) 

Because  we  were  trying  to  represent  discrete  data  (rather  than  continuous  functions),  and 
because  we  were  working  with  regional  rather  than  global  data  sets,  we  felt  that  the  usual 
spherical  harmonic  expansion  techniques  were  not  appropriate.  Instead  we  concentrated  on  the 
use  of  orthogonal  functions  of  discrete  variables. 


We  considered  the  use  of  modified  Chapman  functions  for  representing  altitude  profiles  of 
ion  densities.  These  functions  have  the  advantage  that  peak  height  and  peak  density  are  explicit 
parameters,  but  the  extremely  non-linear  nature  of  these  functions  necessitates  the  use  of  non¬ 
linear  least  squares  fitting  methods.  While  such  methods  produced  excellent  representations  of 
individual  profiles,  the  variation  of  the  fitted  parameters  with  latitude,  longitude  (or  MLT),  and 
UT  was  unacceptably  noisy.  Consequently,  we  chose  to  use  Empirical  Orthomormal  Functions 
for  the  altitude  representation. 

Empirical  Orthonormal  Functions  (EOFs)  have  been  used  extensively  to  represent 
meteorological  and  climatological  data  [Lorenz,  1956;  Kutzbach,  1967;  Davis ,  1976;  and  Peixota 
and  Oort,  1991],  They  have  also  been  used  for  empirical  ionospheric  modeling  [Secan  and 
Tascione,  1984],  EOFs  are  described  in  Appendix  A.  They  have  the  advantage  of  providing  a 
representation  in  terms  of  linear  combinations  of  orthogonal  functions,  which  allows  for 
straightforward  determination  of  coefficients.  However,  because  peak  density  and  peak  height  are 
not  explicit  parameters  of  the  representation,  these  parameters  can  be  determined  only  by 
reconstructing  the  entire  profile  and  invoking  a  peak  finding  algorithm.  We  expect  to  revisit  this 
problem  in  future  versions  of  PIM  and  PRISM  and  implement  a  new  representation  that  combines 
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the  attractive  features  of  both  methods,  i.e.,  that  includes  peak  density  and  peak  height  as  explicit 
parameters  yet  relies  on  linear  combinations  of  orthogonal  functions  to  describe  the  profile  shape. 

For  longitude  (or  local  time)  variations  (and  for  the  low  latitude  F  layer  UT  variation),  the 
obvious  choice  is  a  Fourier  series,  since  trigonometric  functions  retain  their  orthogonality 
properties  on  uniform  discrete  grids  and  because  the  data  is  periodic  in  the  independent  variable. 
These  worked  quite  well  for  the  high  latitude  models  under  all  conditions  and  for  the  low  and 
midlatitude  models  under  low  to  moderate  solar  activity  conditions.  However,  they  did  not  work 
well  for  the  low  and  midlatitude  models  under  high  solar  activity  conditions,  apparently  because 
the  EOF  coefficients  exhibited  exceptionally  large  gradients  at  dawn  and  dusk.  Therefore,  we 
decided  to  tabulate  the  coefficients  in  longitude  for  all  the  low  and  midlatitude  databases. 

For  the  latitude  variations,  we  chose  to  generate  grid-specific  orthogonal  polynomials 
using  the  algorithm  derived  in  Beckmann  [1973]  and  described  in  Appendix  B.  To  help  keep  the 
notation  straight,  we  summarize  it  in  Table  3. 


Table  3.  Notation  Summary 


grid 

variable 

index 

orthogonal  function 

index 

altitude 

z,- 

l<i<  I 

EOF:  S„(z,) 

\<m<M 

latitude 

1  <j<J 

polynomial:  un(X) 

0  <n<N 

longitude 

<P* 

1  <k<K 

trigonometric0:  cos(/7(p)and  sin(ptp) 

0  <p<P 

local  time 

V* 

1  <k<K 

trigonometric*:  cos( p\\f )  and  sin(p\|/) 

0  <p<P 

UT 

1  <1<L 

trigonometric0:  cos(<7x)and  sin(gt) 

0<q<Q 

°Not  used  for  any  database  in  PRISM  1.5,  but  may  be  used  in  future  versions. 
*Used  only  for  high  latitude  databases. 


The  semianalytic  representation  of  each  database  was  generated  in  several  steps.  For  all 
ionospheric  regions,  the  first  step  was  the  determination  of  the  EOFs  from  the  ion  densities  in  the 
database  and  a  set  of  coefficients  clm{Xj,( pj  for  representing  each  ion  density  profile  on  the 
latitude,  longitude,  UT  grid  (See  Appendix  A). 
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M 

n,(  z„X^„,z,)  ..xfrZXl,)  [TDIM] 

m=l 

M 

n,(zl,\l,<ft,ll)  -  [M1DLAT,  ECSD] 

m- 1 
M 

n&Xj, q>*.V,)  -  [LOWLAT] 


where  zn  A,;,  cp*,  t}/*,  t;,  and  \jr,  are  all  points  on  the  model  output  grid,  and  is  the  m* 

EOF  evaluated  at  z;.  (Note,  however,  that  a  different  set  of  gt\z)  functions  are  used  for  each 

ion,  for  each  set  of  geophysical  conditions,  and  for  each  model.) 

For  the  high  latitude  model  (TDIM,  both  E-  and  F-layers),  the  second  step  was  the 
generation  of  Fourier  coefficients  in  MLT,  and  b^{Xj,X for  each  point  on  the 

latitude,  UT  grid. 


(la) 

(lb) 

(lc) 


M  P 


^(zi.,X.,\(/,x/)-£^{^)(X.,x;)co<p\)/)  +  ^)(Xy.,x()sin(p\|r)}gW(z,.)  [TDIM]  (2) 


m= 1  p~ 0 


For  the  low  and  midlatitude  models,  we  found  that  a  truncated  Fourier  series  often  introduced 
spurious  longitudinal  dependences,  apparently  driven  by  the  steep  gradients  at  dawn  and  dusk. 
The  effect  was  particularly  pronounced  at  high  solar  activity  when  the  day/night  contrast  is  the 
greatest.  Consequently,  for  these  models,  the  EOF  coefficients  remain  tabulated  in  longitude. 

For  all  models,  the  next  step  was  the  generation  of  orthogonal  polynomials  from  the 
latitude  grid  (Appendix  B).  For  the  high  latitude  model  (TDIM)  the  coefficients  are  a^(x ,)  and 

P™P(X/),  and  the  ion  density  is  approximated  by 

MNP 

ni(x,.,X,v,x/)  =  ££^{a^(x;)cos(p\}r)  +  (3(^(x/)sin(p\i/)}gm(zi)M„(X)  [TDIM]  (3a) 

m= 1  n=0  p~ 0 


For  MIDLAT  and  ECSD  the  coefficients  arey^X^x,)  and  the  ion  density  is  approximated  by 
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n,(zi,A.,(pt,x,)  =  2X'y^(<P*^/km(z,K(^) 

m~  1  n=0 


[MIDLAT,  ECSD]  (3b) 


For  LOWLAT  the  coefficients  are  7]  ^(<p*,\}L)  and  the  ion  density  is  approximated  by 


U  N 


n,(ZiX  q>*.V/)  “ 


m= 1  n=0 


[LOWLAT]  (3c) 


The  number  of  terms  in  each  series  are  listed  in  Table4  for  each  region. 


Table  4:  Altitude  Grids  and  EOF’s 


Database 

pftnnWssflSflnM 

minimum 

altitude 

maximum 

altitude 

number  of 

EOF’s 

low  latitude  0+ 

55 

160 

1600 

9 

midlatitude  0+ 

49 

125 

1600 

8 

low  &  midlatitude  NO+  &  O2 

28 

90 

400 

7 

high  latitudeO+,  NO+,  &  0*2 

37 

100 

800 

6 

Note  that  in  none  of  these  cases  was  the  altitude  spacing  uniform. 


Because  of  the  extensive  use  of  tabulated  coefficients,  the  ion  density  at  an  arbitrary  point 
must  be  obtained  by  interpolation.  In  PIM  and  PRISM,  altitude  interpolation  is  quadratic,  while 
UT  inteipolation  is  linear.  For  the  MIDLAT  databases,  the  longitude  interpolation  is  also  linear, 
as  is  the  local  time  interpolation  in  the  LOWLAT  databases.  However,  the  longitude  interpolation 
in  the  LOWLAT  databases  is  more  complicated.  First,  the  0+  profile  for  the  desired  magnetic 
latitude  and  local  time  is  reconstructed  for  each  of  the  four  longitude  sectors.  Then  the  peak 
height  and  peak  density  is  determined  for  each  profile.  The  peak  height  for  the  desired  longitude 
is  determined  by  Fourier  interpolation,  and  all  four  profiles  are  shifted  to  match  the  interpolated 
peak  height.  Then  Fourier  interpolation  is  used  at  each  altitude  to  obtain  the  interpolated  ion 
density  profile. 
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3.3  Merging  the  Regional  Models 

Because  we  used  four  different  regional  models  in  the  development  of  PRISM,  the  models 
must  be  merged  at  region  boundaries.  Specifically,  the  low  latitude  and  midlatitude  0+  models 
have  to  be  merged  across  the  boundary  between  low  and  middle  latitudes,  while  all  three  ions 
(0\  NO+,  and  Ot,)  must  be  merged  across  the  boundary  between  midlatitudes  and  high  latitudes. 

The  transition  from  low  latitude  Q+  profiles  to  midlatitude  0+  profiles  takes  place 
between  30°  and  34°  in  both  hemispheres.  The  transition  is  accomplished  by  taking  a  weighted 
average  of  the  hmF2  values  from  the  two  models  in  which  the  weight  shifts  linearly  from  100% 
low  latitude  at  30°  to  100%  midlatitude  at  34°.  The  profiles  are  shifted  to  match  the  averaged 
hmF2  values  and  then  a  similar  weighted  average  of  the  shifted  profiles  is  taken  to  produce  the 
final  merged  profile.  No  transition  for  NQ+  and  O \  is  necessary  since  a  single  model  was  used 
for  these  ions. 

The  transition  from  midlatitude  to  high  latitude  takes  place  over  an  8°  wide  zone  whose 
poleward  boundary  is  the  equatorward  boundary  of  the  trough.  The  transition  process  is  similar 
to  the  low  to  midlatitude  transition,  except  that  the  high  latitude  profiles  are  shifted  to  match  the 
hmF^  and  hmE  values  given  by  the  midlatitude  models.  The  final  profile  is  produced  by  a  weighted 

average  of  midlatitude  and  (shifted)  high  latitude  profiles. 

Although  PIM  and  PRISM  use  geomagnetic  coordinates  internally,  they  can  produce 
output  in  either  geomagnetic  or  geographic  coordinates.  A  contour  map  of  NmF2  in  geographic 
coordinates  (cylindrical  projection)  for  the  June  solstice  at  high  solar  activity  and  moderate 
magnetic  activity  is  displayed  in  Figure  1.  The  equatorial  anomaly  is  clearly  visible  between  East 
longitudes  30°  and  180°,  corresponding  to  local  times  of  1400  and  2400.  The  high  latitude  is 
more  clearly  seen  in  a  polar  projection  such  as  is  displayed  in  Figure  2,  again  in  geographic 
coordinates.  The  figure  shows  the  northern  hemisphere  at  the  December  solstice.  By  is  positive, 

and  the  tongue  of  ionization  resulting  from  a  steady  convection  pattern  is  clearly  visible  on  the 
evening  side. 
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Figure  1.  Contours  of  NmF2  (in  units  of  105  cm'3)  in  cylindrical  equidistant  projection  from  PIM 
for  high  solar  activity,  moderate  magnetic  activity,  at  12  UT  near  the  June  solstice. 
The  equatorial  anomaly  is  clearly  evident  from  about  1400  to  2400  local  time. 
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4.  REAL  TIME  ADJUSTMENT  ALGORITHM 


The  Real  Time  Adjustment  (RTA)  algorithm  for  the  low  and  middle  latitude  regions  is 
different  from  the  algorithm  used  in  the  high  latitude  region.  This  is  partly  due  to  the  relatively 
complex  morphology  of  the  high  latitude  ionosphere  and  partly  due  to  an  evolution  in  our  ideas 
about  the  real  time  adjustment  process  during  the  development  of  PRISM. 

4.1  Available  Data 

The  near  real  time  data  available  for  use  in  the  adjustment  process  comes  from  the  Digital 
Ionospheric  Sounding  System  (DISS),  the  Ionospheric  Monitoring  System  (IMS),  and  a  suite  of 
Special  Sensors  on  the  DMSP  satellites.  The  DISS  network  consists  of  a  network  of  digital 
ionosondes  measuring  critical  frequencies,  critical  heights,  and  bottomside  profiles.  There  are 
projected  to  be  19  DISS  sites  when  the  network  is  complete.  The  Ionospheric  Monitoring 
System  consists  of  a  separate  network  of  dual  frequency  GPS  receivers  measuring  TEC.  The 
number  is  IMS  sites  is  not  finallized  but  is  expected  to  be  about  five.  The  DMSP  Special  Sensors 
include  the  in  situ  plasma  properties  measured  by  the  SSIES  instrument  and  the  precipitating 
particle  measurements  of  the  SSJ/4  instrument.  In  the  future,  the  SSJ/4  instrument  will  be 
replaced  by  the  SSJ/5  instrument  and  a  new  set  of  ionospheric  remote  sensing  instruments,  SSUSI 
and  SSULI  will  be  added.  SSUSI  (a  multispectral  UV  imager)  and  SSULI  (a  multispectral  limb 
imager)  will  make  measurements  of  dayglow,  nightglow,  and  auroral  optical  emissions  (mostly 
ultraviolet).  The  observed  intensities  will  be  processed  on  the  ground  to  deduce  ionospheric 
properties  in  the  form  of  PRISM  profile  parameters  [e.g..  Fox  et  al.,  1994],  The  two  UV  imagers 
will  be  flown  on  DMSP  beginning  some  time  near  the  end  of  the  decade. 

Because  TEC  is  an  integral  quantity,  it  is  not  easy  to  incorporate  it  into  the  PRISM  real 
time  scheme.  For  PRISM  1.5,  we  have  chosen  to  convert  TEC  into  an  equivalent  point  datum. 
Our  method  for  doing  so  is  described  in  Appendix  C. 

4.2  Low  and  Midlatitude  Adjustment  Parameters 

The  PRISM  real  time  adjustment  algorithm  operates  on  six  parameters  that  prescribe  how 
an  electron  density  profile  is  to  be  modified  or  "corrected": 
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1.  A f0F2 ,  the  correction  to  the  model  faF2, 

2.  A faE,  the  correction  to  the  model  faE, 


3.  A hmF2 ,  the  correction  to  the  model  hmF2 , 

4.  A hmE,  the  correction  to  the  model  hmE, 

5.  AN,op ,  the  correction  to  the  0+  density  at  a  specific  altitude  (i.e.,  the  DMSP  altitude) 


6.  A Htop,  the  correction  to  the  0+  scale  height  at  a  specific  altitude  (i.e.,  the  DMSP 
altitude) 


The  nominal  value  for  each  of  these  parameters  is  zero.  A  positive  (negative)  value  means  that 
the  model  value  must  be  increased  (decreased).  Using  the  available  near  real  time  data,  the  real 
time  adjustment  process  will  assign  non- zero  values  at  each  location  where  data  is  available  (the 
driver  sites). 

Parameters  1-4  are  based  on  direct  measurements  by  DISS  digital  ionosondes.  Parameter 
5  is  based  on  the  direct  measurement  of  the  0+  density  by  the  SSIES  instrument  on  board  the 
DMSP  satellite.  In  contrast,  Parameter  6  must  be  inferred  from  the  electron  and  ion  temperatures 
(Te,  Ti )  measured  by  the  SSIES  instrument.  At  midlatitudes,  PRISM  assumes  that  the  topside 
ionosphere  is  in  diffusive  equilibrium  and  calculates  the  topside  scale  height  from 


k(Te+T,) 


m0.g 


(5) 


The  corresponding  model  value  is  obtained  from  the  model  0+  densities  (n2,  nt)  at  the  grid 
altitudes  (z2,  z})  immediately  above  and  below  the  DMSP  altitude: 

(6) 


In  {njn2) 
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At  low  latitudes,  the  topside  scale  height  is  determined  by  diffusive  processes  alone,  so  the  scale 
height  cannot  be  inferred  from  the  SSIES  temperature  data. 


4.3  Real  Time  Adjustment  of  the  Low  and  Midlatitude  Profile  Parameters 

In  PRISM  1.5,  after  the  values  of  the  profile  adjustment  parameters  have  been  determined 
at  each  driver  site,  the  global  correction  field  is  determined  using  a  weighted  average  method. 

1.  Given  N  driver  locations  and  the  associated  geomagnetic  coordinates  (Ln,§n),  let  the 
unadjusted  PRISM  value  at  the  nth  point  be  un  and  the  measured  value  be  v„. 

2.  The  correction  to  be  applied  to  the  unadjusted  PRISM  (i.e.,  PIM)  value  at  the  n,h 
point  is  cn=v„-un. 

3.  At  any  other  point,  (^,<J)),  the  correction  to  be  applied  to  the  PIM  value  is 

c(M>)=-^ -  (?) 

n- 1 

where  the  are  weight  functions  that  depend  on  the  distance  measure  dn(X,ty) 

between  the  point  and  the  n'h  point(A, „,<()„). 

■-.(>■■♦)=  ™  e  (8) 

dn{L,Q>)  ££ 

In  PRISM  1.5  we  have  used  the  following  distance  measure. 

-cosy  „(^-.<l>)]  (9) 

where  Y„(X,<|>)  is  the  great  circle  distance  between  (Xd)  and  (?Ln,())J: 
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cos y  „  ( A. ,  <|) )  =  sin  X  sin  A,  n  +  cos  A,  cos  A.  n  cos(<()  -  <|)  „ ) 


(10) 


As  a  practical  matter,  the  actual  weighting  function  used  in  PRISM  is  wa(A,,<p)  =  \/dn{X,g>) 
unless  dm(A,,cp)  <8  for  some  driver  station  m.  In  that  case,  wm  =  land  wn  =  0  for  n  -f-  m.  In 
PRISM  1.5,  8  is  set  to  1CT*. 

This  method  ensures  that  PRISM  will  reproduce  the  input  data.  However,  because  the 
decorrelation  length  of  the  ionosphere  is  of  the  order  of  1000  km  or  less,  no  interpolation  scheme 
can  hope  to  accurately  reproduce  ionospheric  parameters  where  there  is  no  data.  Clearly,  the 
denser  the  data  net,  the  better  the  model  will  do.  Unfortunately,  at  least  initially,  the  data  will  be 
quite  sparse.  In  order  to  ensure  that  sparse  datasets  do  not  produce  unreasonable  ionospheric 
specifications,  PRISM  1.5  checks  the  spatial  distribution  of  the  data  and  inserts  “phantom 
stations”  with  the  corrections  forced  to  zero.  The  procedure  for  placing  phantom  stations  is  as 
follows. 

The  low  and  midlatitude  region  is  divided  into  32  rectangular  subregions  with  boundaries 
defined  by  the  following  parallels  of  latitude  and  meridians  of  longitude: 

parallels  of  latitude:  60°S,  30 °S,  0°,  30°N,  and  60°N 

meridians  of  longitude:  30°E,  75°E,  120°E,  165°E,  210°E,  255°E,  300°E,  and 

345°E 

There  is  a  phantom  station  located  at  the  center  of  each  rectangle.  During  the  data  ingestion 
process,  PRISM  keeps  track  of  how  many  stations  are  located  in  each  rectangular  subregion,  and 
how  many  have  f0F2  data.  As  long  as  at  least  one  station  in  a  subregion  reports  a  value  for  f0F2, 
then  the  phantom  station  for  that  subregion  is  ignored.  If,  however,  no  station  in  the  subregion 
reports  a  value  for  fQF2,  then  the  phantom  station  is  assigned  a  value  equal  to  the  PIM  value  at 
that  point.  A  similar  procedure  is  followed  for  hmF2,  FaE,  and  hmE  data.  These  phantom  stations 
force  the  correction  field  to  relax  to  values  near  zero  in  regions  where  there  is  no  data. 

4.4  Modifying  the  Low  and  Midlatitude  Model  Profiles 

The  interpolation  method  described  above  provides  a  global  correction  field  that  can  be 
used  to  calculate  the  profile  correction  parameters  at  any  location.  In  this  section,  we  describe  the 
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way  in  which  the  profiles  are  adjusted  using  the  profile  correction  parameters.  First  the  layer 
heights  are  adjusted,  then  the  layer  peak  densities,  and  finally  the  topside  correction  is  applied. 

The  layer  height  correction,  using  parameters  3  and  4  is  a  simple  shifting  of  the  profiles.  If 
the  F  layer  correction  is  A hmF2 ,  then  the  0+  profile  is  shifted  so  that 

%r(z)  =  n0.°“(z-Ah.F2)  (11) 

A  similar  shift  is  applied  to  the  molecular  ion  profiles,  except  that  the  altitude  shift  is ^hmE. 

Once  the  altitude  corrections  have  been  applied,  parameters  1  and  2,  the  critical 
frequencies,  will  be  used  to  scale  the  ion  density  profiles.  For  the  F  layer  (0+),  the  unadjusted 
peak  density  is  converted  to  a  critical  frequency.  The  correction  to  be  applied  (determined  by  the 
algorithm  in  section  4.3)  is  simply  added  to  the  unadjusted  frequency.  The  corrected  frequency  is 
then  converted  back  to  a  density.  The  ratio  of  the  corrected  density  to  the  unadjusted  density  is 
used  to  scale  the  profile  above  the  peak.  Below  the  peak,  an  additive  correction  is  applied.  This 
additive  correction  smoothly  vanishes  as  the  altitude  approaches  hmE.  This  ensures  that  F  layer 
corrections  do  not  impact  the  E  layer.  This  additive  correction  is  described  in  Appendix  C.  A 
similar  scaling  is  applied  to  the  molecular  ions  except  that  the  additive  correction  is  used  above 
the  peak  and  smoothly  vanishes  as  the  altitude  approaches  hmF2.  This  ensures  that  the  E  layer 
corrections  do  not  impact  the  F  layer. 

Parameters  5  and  6  are  used  to  correct  the  topside  profile  based  on  ne,  n{,  Te,  and  T{ 
measurements  from  SSIES  on  DMSP  (nominally  840  km).  Let  Nm(z )  be  the  model  0+  profile 
(after  fj7 2  and  hnf2  corrections  have  been  applied).  Further,  let  zp  =hmF2,  and  let  ztop=  the 

altitude  at  which  AN,op  and  A H,op  were  measured  (usually  the  DMSP  altitude).  Let  Nc(z)  be  the 
corrected  profile  based  on  A Ntop  and  A Htop.  Nc(z )  must  satisfy  the  following  conditions: 


Nc{",op)  =  N,op  =  Nm(z,op)  +  ANtop 


(12) 


and 


1  dNc 
AL  dz 


L  Htop  Hm(zJ  +  Mitop 


(13) 
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In  PRISM  the  corrected  profile  is  obtained  from  the  model  profile  by  scaling  the  altitude  so  that 


where 


1 
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(14) 


(15) 

(16) 


(17) 


Direct  substitution  of  z  =  ztop  in  Equation  (15)  verifies  that  C,(ztop)  =  z1  so  that  Equation  (12)  is 
satisfied.  That  Equation  (13)  is  satisfied  may  be  seen  from 
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1  dN. 
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(18) 


This  form  ensures  that  the  correction  and  its  first  derivative  vanish  at  the  peak  and  that  the  0+ 
concentration  and  its  slope  are  as  specified  at  z,  . 

When  there  is  no  topside  scale  height  data  (as  at  low  latitudes),  PRISM  sets  Hlop  to  the 

value 

z,  —  z 

H  = - --//m  (when  Ht  is  not  determined  from  data)  (19) 

z  —  Z 

top  ^p 
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so  that  the  altitude  scaling  becomes  linear 


C  (z)  =  z,+ 


(when  Htop  is  not  determined  from  data) 


(20) 


An  example  of  the  midlatitude  profile  correction  algorithm  applied  to  a  single  profile  is 
shown  in  Figure  3,  for  which  we  have  used  Incoherent  Scatter  Radar  (ISR)  data  from  Arecibo  to 
simulate  simultaneous  digisonde  and  SSIES  data.  In  this  case,  the  data  is  from  1800  UT  on  4 
October  1989.  In  Figure  3a,  the  ISR  data  is  compared  with  the  climatology  of  PIM  (i.e., 
unadjusted  PRISM).  A  profile  from  the  Ionospheric  Conductivity  and  Electron  Density  (ICED) 
model  [T as  done  et  al.,  1988],  which  is  currently  operational  at  AFSFC  is  also  shown.  The  actual 
NmF2  is  a  factor  of  two  higher,  and  the  actual  hmF2  is  about  40  km  higher,  than  the  climatological 
values.  In  Figure  3b,  the  PRISM  profile  has  been  adjusted  only  to  match  the  actual  NmF2.  In 
Figure  3c,  both  NmF2  and  hmF2  have  been  adjusted  to  match  the  actual  data.  Finally,  in  Figure  3d, 
both  the  peak  and  topside  adjustments  have  been  made.  If  we  define  an  RMS  density  error  as 


RMS  density  error  = 


2 

1  13 
— y 

N  PRlSAf(Zk  )  ~  ^  data  ( Zk  ) 

i  o  ^ 
k—\ 

Ndata{Zk) 

(21) 


then  PRISM’s  RMS  density  error  declines  from  46%  with  no  profile  adjustment  to  25%  with  an 
NmF2  adjustment  only,  to  18%  with  both  NmF2  and  hmF2  adjusted,  and  to  16%  with  both  peak  and 
topside  adjustments. 

Note  that  in  operational  use  there  are  normally  many  data  sites,  and  the  topside  and 
bottomside  data  are  seldom  colocated.  Operationally,  PRISM  uses  all  available  hmF2  and  hmE 
data  to  establish  global  layer  height  correction  fields,  then  it  uses  the  available  NmF2  and  NmE 
data  to  calculate  global  profile  scaling  fields,  and  finally  it  uses  the  topside  data  to  calculate  global 
topside  correction  fields.  These  correction  fields  are  then  applied  to  each  profile  on  the  output 
grid  as  described  above. 
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Figure  3.  The  PRISM  profile  adjustment  procedure  illustrated  using  Incoherent  Scatter  Radar 
(ISR)  data  from  Arecibo.  (a)  No  adjustment  (i.e.,  PIM  profile):  The  RMS  density 
error  is  46%.  (b)  NmF2  adjustment  only:  The  RMS  density  error  is  25%.  (c)  NmF2 
and  hmF2  adjustment:  The  RMS  density  error  is  18%.  (d)  NmF2,  hmF2,  and  topside 
adjustments:  The  RMS  density  error  is  16%. 
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4.5  The  High  Latitude  Adjustment  Algorithm 

Due  to  the  complexity  of  the  high  latitude  ionosphere,  the  real  time  adjustment  algorithm 
differs  appreciably  from  the  low  and  midlatitude  algorithm.  Until  SSUSI  data  becomes  available, 
there  will  be  insufficient  data  to  adjust  the  parameterized  USU  model  in  the  way  that  the 
parameterized  low  and  midlatitude  models  can  be  adjusted.  Even  with  SSUSI  data,  it  is  not  clear 
that  the  midlatitude  algorithm  is  appropriate  for  the  high  latitude  regions. 

In  PRISM  1.2,  the  first  step  in  the  high  latitude  real  time  adjustment  process  is  the 
establishment  of  boundary  locations.  Three  boundaries  are  required:  (1)  the  equatorward  edge  of 
the  trough,  (2)  the  equatorward  edge  of  the  auroral  oval,  and  (3)  the  poleward  edge  of  the  auroral 
oval. 

The  equatorward  edge  of  the  trough  is  determined  from  SSIES  drift  meter  data  as  the 
point  where  the  measured  ion  drift  speed  departs  from  the  corotation  value.  The  trough  boundary 
as  a  function  of  magnetic  longitude  is  given  by  the  formula 


e,((p)=60(/,)  +  aexp 


r<p-6Y 

V  C  J 


where  (p  is  the  magnetic  local  time  (MLT,  hours),  and  7t  is  a  "trough  index"  correlated  with  Kp. 
0O(/,)  is  the  radius  of  the  trough  boundary  at  magnetic  local  midnight  and  is  given  by 

0o(/,)  =  24.4°  +  2.12°/, 


The  second  term  represents  the  dayside  distortion  of  the  boundary,  which  would  otherwise  be  a 
circle  centered  on  the  magnetic  pole.  The  parameter  values  are 

a  =  -10.5° 
b  =  11.5  hr 
c  =  3.88  hr 
p  =  2.73 
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This  model  is  an  approximation  to  the  convection  boundaries  shown  in  Heppner  and  Maynard 
[1987]. 

The  trough  index,  /,,  is  determined  according  to  the  following  algorithm. 

1.  If  there  is  no  ion  drift  data,  or  if  a  boundary  cannot  be  identified  in  the  data,  then/,  =  Kp. 

2.  If  a  single  boundary  crossing  is  identified  at  colatitude06  and  local  timetp^ ,  then 

0O  =0, -aexp 
and 

j  e0-24.4° 

'  2.12° 

3.  If  two  or  more  crossings  are  identified,  then  the  value  of  /,  used  is  the  average  of  the  values 
determined  for  each  crossing. 

The  boundaries  of  the  auroral  oval  are  determined  from  electron  and  ion  precipitation  data 
from  the  SSJ/4  instrument.  Separate  boundaries  are  determined  for  electron  and  ion  precipitation. 
The  algorithm  for  determining  the  electron  and  ion  precipitation  boundaries  from  the  SSJ/4  data  is 
described  in  Appendix  D. 

The  next  step  in  the  high  latitude  adjustment  process  depends  on  the  amount  and  kind  of 
data  available.  The  decision  of  how  to  proceed  is  made  separately  for  the  E  layer  (NO+and  Oj) 
and  the  F  layer  (0+).  In  each  case,  two  choices  are  available: 

F  layer: 

1.  Perform  a  simple  least  squares  adjustment  of  theUSU  0+  model. 

2.  Use  a  semi-empirical//^  model  (FMQDEL)  to  adjust  the  USU  0+  profiles. 

E  layer: 
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1.  Perform  a  simple  least  squares  adjustment  of  the  USUNO+  and  0*2  models. 

2.  Use  a  fast,  first  principles,  E  layer  local  chemistry  model  (HLE). 

The  decision  matrix  used  in  PRISM  is  shown  in  Table5. 


Table  5.  PRISM  High  Latitude  Decision  Matrix 


|  Data  Available? 

Model  Used  by  PRISM  | 

DISS 

Yes 

Yes 

Yes 

FMODEL 

HLE 

Yes 

Yes 

No 

FMODEL 

USU 

Yes 

No 

Yes 

FMODEL 

HLE 

Yes 

No 

No 

USU 

USU 

No 

Yes 

Yes 

USU 

USU 

No 

Yes 

No 

USU 

USU 

No 

No 

Yes 

USU 

USU 

No 

No 

No 

USU 

USU 

In  experimenting  with  high  latitude  data,  we  found  that  extrapolating  SSJ/4  data  taken 
along  the  DMSP  orbital  track  to  points  well  away  from  the  orbital  track  was  very  risky.  We 
found  no  suitable  model  of  the  instantaneous  auroral  precipitation  for  this  extrapolation.  When 
SSUSI  auroral  image  data  becomes  available,  this  limitation  will  be  removed  because  much  less 
extrapolation  will  be  required.  It  should  be  possible  to  use  HLE  whenever  timely  SSUSI  images 
are  available. 

FMODEL  is  a  semi-empirical  model  of  faF2  based  on  a  combination  of  theory  and  data.  It 
is  divided  into  three  regions:  the  subauroral  trough,  the  auroral  oval,  and  the  polar  cap.  The  f0F2 

determined  by  least  squares  adjustment  of  the  model  parameters  is  used  to  scale  the  USU  0+ 
profiles.  No  further  adjustment  of  the  profiles  is  performed. 

The  subauroral  trough  is  divided  into  two  local  time  regimes:  evening  (from  1200  to  0000 
MLT)  and  morning  (from  0000  to  1200  MLT).  Each  trough  (morning  and  evening)  has  a  depth 
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parameter  that  specifies  the  difference  between  the  midlatitude  value  of  faF2  at  the  equatorward 
edge  and  the  trough  minimum.  The  local  time  variation  of  the  trough  depth  is  fixed  (not  part  of 
the  least  squares  adjustment  process).  If  the  width  of  the  trough  is  less  than  3°,  the  depth  is 
reduced  in  proportion  to  the  width  so  that  when  the  width  vanishes  so  does  the  depth.  The 
thickness  of  the  (equatorward)  trough  wall  is  always  60%  of  the  total  width  of  the  trough.  The 
poleward  edge  of  the  trough  is  the  equatorward  edge  of  the  auroral  F-layer,  so  the  poleward 
"wall"  is  considered  to  be  part  of  the  auroral  region. 

At  fixed  magnetic  local  time,  the  auroral  F  layer  fJF2  is  simply  a  cubic  polynomial  in 
magnetic  latitude: 


fA  W  =  /_  +  +B(X  -  X_J 

where  fmax  =  f0F2('^m/,x)  is  an  extremum,  and  A  and  B  are  chosen  so  that  faF2  is  continuous 
across  the  boundaries  with  the  trough  and  polar  cap. 

The  background  polar  cap  faF2  is  obtained  from  the  URSI  coefficients  using  an  internally 
derived  effective  sunspot  number,  the  value  of  which  is  determined  as  part  of  the  least  squares 
adjustment  process  based  on  iono  sonde  data.  It  should  be  noted  that  this  part  of  the  model 
describes  only  the  background  polar  cap  ionosphere,  not  "polar  cap  patches”  or  "polar  cap  arcs." 
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5.  VALIDATION 


Prior  to  delivery  to  the  Air  Force  Space  Forecast  Center,  and  earlier  version  of  PRISM 
(1.2)  underwent  an  extensive  validation  using  historical  data  approximating  the  kinds  of  data  that 
will  be  available  to  PRISM  in  operational  use.  The  results  of  this  validation  effort  were  reported 
in  Daniell  et  al.  [1994]  and  are  not  repeated  here.  Since  PRISM  1.5  incorporates  several 
substantial  changes  from  PRISM  1.2,  a  new  validation  will  be  carried  out  in  the  future  under  a 
separate  contract. 
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6.  DISCUSSION 


We  have  described  the  development  of  PRISM,  a  real  time  ionospheric  specification 
model  based  on  parameterized  theoretical  ionospheric  models.  Unlike  previous  specification 
models,  PRISM  is  based  on  physical  models  rather  than  empirical  models.  PRISM  consists  of 
two  parts:  the  parameterized  physical  models  (PIM)  and  a  real  time  update  system  that  ingests 
both  ground-based  and  space-based  data  and  modifies  PIM  profiles  accordingly.  PRISM  provides 
considerable  improvement  over  simple  climatology  in  the  vicinity  of  the  data  sources  and  does  no 
worse  than  climatology  at  locations  remote  from  the  data.  With  the  advent  of  remote  ionospheric 
sensing  using  UV  airglow  and  auroral  emissions  measured  from  satellites  and  the  proliferation  of 
dual  frequency  GPS  receivers  providing  line  of  sight  TEC  measurements,  we  expect  this  approach 
to  provide  improved  ionospheric  specifications  in  near  real  time  for  the  AFSFC  and  its  customers. 

A  number  of  compromises  were  required  in  the  development  of  PIM  and  PRISM.  First, 
of  course,  was  the  necessity  of  using  parameterizations  (in  the  form  of  diumally  reproducible 
runs)  of  the  physical  models,  rather  than  the  physical  models  themselves.  Second,  we  had  to  use 
empirical  models  (e.g.,  MSIS-86)  instead  of  physical  models  to  provide  the  necessary  inputs  to 
the  ionospheric  models.  Third,  we  had  to  use  a  tilted  dipole  representation  of  the  earth’s 
magnetic  field  instead  of  a  more  realistic  model.  This  last  compromise  is  mitigated  somewhat  by 
the  use  of  Corrected  Geomagnetic  (CGM)  coordinates  in  PIM  and  PRISM.  While  not  fully  self- 
consistent,  this  does  allow  a  more  realistic  representation  of  geomagnetically  controlled  features 
such  as  the  equatorial  anomaly.  We  expect  to  remove  the  compromises  as  available  computing 
power  increases  in  the  future. 

The  particular  features  described  here  apply  to  PRISM  version  1.5  and  PIM  version  1.3. 
Significant  enhancements  to  both  PIM  and  PRISM  are  planned  for  the  near  future.  H+  ion 
densities  based  on  a  parameterization  of  the  plasmasphere  model  of  Bailey  and  Sellek  [1990]  will 
be  added  so  that  PIM  and  PRISM  can  give  electron  density  profiles  up  to  the  plasmapause.  At 
the  same  time  the  coefficient  files  will  be  regenerated  using  a  single  global  ionospheric  model, 
eliminating  the  need  to  merge  models  across  region  boundaries.  The  parameterization  process 
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will  also  be  reexamined  in  order  to  produce  a  more  accurate  and  more  efficient  analytic  fit  to  the 
model  runs.  The  resulting  version  of  PRISM,  while  considerably  enhanced,  will  be  designed  to  fit 
into  AFSFC’s  operational  configuration  with  no  changes  to  the  existing  software  that  interacts 
with  PRISM  1.5. 
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Appendix  A,  Empirical!  Orthonormal  Functions 

This  treatment  of  empirical  orthogonal  functions  (EOF’s)  is  based  on  the  Appendix  of 
Secan  and  Tascione  [1984],  which  was  based  on  Lorenz  [1956],  Kutzbach  [1967],  and  Davis 
[1976].  See  also  Peixota  and  Oort  [1991].  The  reader  is  referred  to  these  references  for 
mathematical  proofs  of  the  assertions  made  below.  In  the  following  discussion,  we  use  the 
notation  given  in  Table  3  of  the  main  text. 

A  database  consists  of  altitude  profiles  at  certain  longitudes,  certain  latitudes,  and  certain 
Universal  Times.  (See  Tables  1-4  of  the  main  text.)  Let  S  be  the  number  of  altitude  profiles  in  a 
database,  and  let  I  be  the  number  of  points  in  each  altitude  profile.  We  would  like  to  represent 
each  altitude  profile  of  the  quantity  "F  (e.g.,  0+ concentration)  as  an  expansion  in  orthogonal 
functions,  gm(z;): 


M 


^•)  =  Xa^i)  +  ^U).  s  =  l...S,  i  -  1...I 


(Al) 


m= 1 


where  rs(z;)  is  the  residual,  and  the  coefficients®^  are  calculated  from 

(A2) 

1=1 

In  principle,  any  orthogonal  set  of  functions  may  be  used.  However,  the  references  cited  above 
provide  an  algorithm  for  finding  the  set  which  minimizes  the  RMS  error  for  a  given  number  of 
terms,  M  <  /.  We  summarize  the  algorithm  here. 

First  define  the  /  x  /  covariance  matrix  C  with  elements 

(Al) 

^  5=1 

Now  consider  the  eigenvalue/eigenvector  problem C(p  =  <pL  or 
/  / 

EQP  *  =  5X5  A  =  <P  A  (A2) 

;'= i  J= i 

where  (p  =  {(pyj  is  the  matrix  of  eigenvectors  of  C  =  [Cy- j,  and  L  =  [5yA,y]  is  a  diagonal  matrix 
whose  elements  are  the  corresponding  eigenvalues.  (The  column  of  (p  is  the  eigenvector 
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corresponding  to  the  eigenvalue,  Xk.)  By  convention,  the  eigenvectors  and  eigenvalues  are 
ordered  so  that  X1>X2>...>XI.  Because  C  is  a  real  symmetric  matrix,  eigenvectors 
corresponding  to  unique  eigenvalues  are  guaranteed  to  be  orthogonal  [See,  e.g.,  Hildebrand , 
1965].  Because  of  the  origin  of  the  matrix  C,  it  is  unlikely  that  any  of  its  eigenvalues  will  be 
degenerate,  so  we  may  assume  that  (p  is  an  orthogonal  set.  According  to  Secan  and  Tascione 
[1984]  and  references  therein,  the  set  of  orthogonal  functions  that  minimizes  the  RMS  error  for  M 
terms  is  just  the  first  M  eigenvectors: 

i'  =  l,2,...7;  m  =  \,2,...M  (A3) 

These  are  the  Empirical  Orthonormal  Functions  (EOFs). 

As  a  practical  matter,  we  have  found  that  the  number  of  EOFs  needed  to  provide  a 
reasonably  good  representation  for  all  the  profiles  is  about  7/6,  as  illustrated  in  Table  4  in  the  main 
text.  The  only  exception  is  the  low  and  midlatitude  E  layer  (NO+  and  0\),  probably  because 
these  databases  covered  both  hemispheres  simultaneously.  We  have  also  found  that  substantial 
improvement  in  representation  does  not  occur  until  the  number  of  EOF’s  is  about  7/2. 
Furthermore,  the  EOF’s  derived  for  one  database  were  inadequate  for  any  other  database,  and  the 
EOF’s  simultaneously  derived  from  several  databases  produce  noticeably  poorer  representations 
than  those  derived  for  each  database  individually.  Consequently,  we  have  derived  separate  EOF 
sets  for  each  database. 

The  first  nine  EOF’s  derived  from  the  low  latitude  F  region  (0+)  database  for  the  US 
longitude  sector,  the  December  solstice,  and  moderate  solar  activity,  are  shown  in  Figure  Al. 
The  first  EOF  always  has  the  least  structure,  and  successive  EOF’s  become  progressively  more 
structured.  Although  differing  in  detail,  the  EOF’s  for  the  other  databases  are  qualitatively 
similar. 
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Altitude  (km) 


1800 


Low  Latitude  0 


Figure  Al.  The  Empirical  Orthonormal  Functions  (EOF’s)  for  low  latitude  0+  derived  from  the 
LOWLAT  output  databases  for  the  USA  longitude  sector,  December  solstice, 
moderate  magnetic  activity,  and  moderate  solar  activity.  Only  the  first  nine  EOF’s 
are  plotted  because  these  are  the  ones  used  in  PIM  and  PRISM. 
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Appendix  B.  Orthogonal  Polynomials  of  Discrete  Variables 


Because  the  databases  to  which  we  desire  analytic  approximations  have  discrete  latitude 
grids,  we  preferred  to  use  polynomials  whose  orthogonality  is  defined  in  terms  of  that  grid,  rather 
than  in  terms  of  integrals  over  the  interval.  The  algorithm  for  generating  orthogonal  polynomials 
on  a  specified  grid  is  given  by  Beckmann  [1973].  Let  us  denote  the  desired  polynomials  by  un(X) 

and  define  u. ,(A.)  =  0  and  u0(X)  =  1.  Note  that  the  polynomials  are  continuous  functions  of  the 
continuous  variable  X  even  though  their  orthogonality  is  defined  in  terms  of  the  discrete  grid 
[ Xj,j  =  1,2,... /j.  The  recursion  relation  for  the  polynomials  is 

<bi) 

“n-l 

where  the  norms  hn  are  given  by 

*.=5>.!(7-,)  ®2> 

and  the  recursion  constants  Bn  are  given  by 

B .  =irXh“.N  <B3) 

Hn  7=1 

The  polynomials  generated  by  this  algorithm  may  be  used  to  represent  the  latitude 
variations  of  the  Fourier  coefficients  amp  and  bmp  (see  main  text): 

amp{^j^l)  =  Y<a'nnp{T‘l)Un(^i)  (B4) 

n~l 

l)  =  l)Un(^j)  (D5) 

n=l 

where 

fimnpfo  l)  =TT^iamp(^j'>X  l^n^j) 

nn  7=1 

nn  7=1 
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Appendix  C.  F-8ayer  Density  Scaling  Algorithm  and  TEC  Adjustments 

In  order  to  be  able  to  adjust  E-  and  F-layer  parameters  separately,  PRISM  uses  an  altitude 
dependent  scaling  algorithm  for  each  layer.  Because  the  TEC  adjustment  process  is  dependent  on 
the  details  of  the  F-layer  adjustment  process,  they  are  described  together. 

For  altitudes  above  hmF2  the  adjustment  is  a  simple  multiplicative  scaling: 


where 


(old) 


z  >  hR 


(Cl) 


SF  — 


NF2{old) 


(C2) 


When  is  derived  from  DISS  or  other  ionosonde  data,  it  is  simply 


NmF±new)  =  1.24x!04(/oF2(diss))2 


(C3) 


where  faF2D]SS)  is  the  measured  value  in  MHz  and NmF2nen,)  is  incm'3. 

For  altitudes  below  hmF2,  the  adjustment  is  additive  and  vanishes  gradually  as  z 
approaches  hmE. 


n£"(z)  =  n™(z)+An0.(z) 

where 


(C4) 


A nQ.(z)  = 


exp|  2  l~exPi 


hmE  —  hmF2 

m _ m  l 

W 


.  (hmE-hR 

1-exp  — - ^ 

*1  W 


0, 


(C5) 


z  <  hE 
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The  value  of  the  parameter  W  is  determined  by  the  shape  of  the  bottomside  0+  density  profile, 
and  is  usually  of  the  order  of  100  km. 

In  order  to  use  TEC  data,  PRISM  must  convert  the  TEC  measurement  into  an  equivalent 
point  measurement.  TEC  data  is  ingested  as  vertical  equivalent  TEC  at  the  Ionospheric 
Intersection  Point  (IIP).  If  the  IIP  is  farther  than  1000  km  from  the  nearest  DISS  site,  then  the 
TEC  is  converted  into  an  equivalent  NmF2.  PRISM  converts  the  observed  vertical  equivalent 
TEC  (TECobs)  into  an  F-layer  correction  factor  that  is  applied  in  the  same  way  as  the  correction 
factors  derived  from  DISS  (ionosonde)  data.  First,  PRISM  calculates  a  TEC  correction  as 

ATEC  =  TECobs  —  TECprism  (C6) 

The  equivalent  NmF2  correction  is  determined  from 


or 


where 


ATEC=fK”')(z)-'!o“)w]* 

ATEC  =  WbANmF2  +SFj^n^\z)d2 

Wt 


Therefore,  the  value  of  NmF^new'>  corresponding  to  TECobs  is 


TEC,te-TECMSM-iytA'„FrlA,  w 


WbNmFf 


JK,Fi  n°* 


(C7) 

(C8) 

(C9) 

(CIO) 


From  this  point,  PRISM  treats  TEC  data  as  if  it  were  ionosonde  data  with  hmF2  missing. 


If  the  IIP  is  within  1000  km  of  a  DISS  site,  then  the  TEC  measurement  is  converted  into 
an  equivalent  SSIES  density  measurement.  First,  the  Af0F2  and  A hmF2  corrections  determined 
from  the  nearest  DISS  site  is  applied  to  the  PIM  profile  at  the  IIP.  Then  the  topside  portion  of 
the  profile  is  corrected  to  force  agreement  with  the  TEC  measurement: 

(cii> 


(z,  z  <  hmF2 

\hmF2  +a(z  -  hmF2),  z  >  hmF2 


(Cl  2) 


where  a  is  a  scale  factor  to  be  determined  from  TEC. 


The  difference  in  TEC,  as  defined  by  Equation  fi),  is 

ATEC  =  f“  n("r]dz-  f  n(°u)dz 

JhmF2  O  O 


(Cl  3) 


or 


ATEC  = 


resulting  in 


a  = 


(Cl  4) 


(Cl  5) 


If  we  identify  the  scale  factor#  with  the  ratio  (z;  ~hmF2)J(zd  ~hmF2),thcn 


zi=KF2+{zd-KF2)a 


(Cl  6) 
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N  top  =  «o°U)(Zl) 


(Cl  7) 


and 

H,op=aHm  (CIS) 

If  the  IIP  is  within  1000  km  of  both  a  DISS  station  and  a  DMSP  orbital  track,  the  TEC 
measurement  is  ignored.  In  a  future  version  of  PRISM,  in  which  the  plasmasphere  is  included, 
TEC  measurements  in  such  situations  will  be  used  to  constrain  theH+  density. 
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Appendix  B.  Auroral  boundary  determination 

There  are  four  auroral  precipitation  boundaries  in  each  hemisphere:  The  equatorward  and 
poleward  electron  boundaries  and  the  equatorward  and  poleward  ion  boundaries.  Since  there  are 
always  two  DMSP  satellites  in  orbit,  there  are  potentially  four  crossings  of  each  boundary.  The 
specification  of  a  given  boundary  depends  on  the  number  of  crossings  detected  in  the  data.  The 
boundary  detection  is  performed  as  follows. 

1.  If  SSJ/4  data  are.  available,  then  divide  each  orbit’s  worth  of  data  into  four  segments,  each 
segment  extending  from  the  most  equatorward  point  to  the  most  poleward  point 

2.  Each  segment  is  searched  from  equator  toward  the  pole  until  the  electron  energy  flux  exceeds 
0.25  erg  cm'2  s'1  (for  the  equatorward  electron  boundary)  and  the  ion  energy  flux  exceeds 
0.1  erg  cm'2  s'1  (for  the  equatorward  ion  boundary). 

3.  Search  each  segment  from  pole  toward  the  equator  until  the  same  thresholds  are  exceeded  to 
establish  the  poleward  boundaries. 

The  boundaries  are  always  assumed  to  be  circular.  The  specification  of  each  boundary  circle  (in 
terms  of  its  center  and  radius)  depends  on  the  number  of  crossings  detected. 

1.  If  four  crossings  are  found,  there  are  four  possible  combinations  of  three  points.  Each  of  the 
four  possible  combinations  defines  a  circle,  which  may  be  described  by  its  center  and  radius. 
The  center  coordinates  and  radii  of  the  four  circles  are  averaged  to  determine  the  boundary 
used  in  PRISM. 

5.  If  only  three  crossings  are  found,  the  boundary  is  uniquely  defined  by  the  circle  passing 
through  the  three  points. 

6.  If  only  two  crossings  are  found,  the  boundary  is  taken  to  be  the  circle  that  passes  through  both 
points  and  has  a  radius  equal  to  the  mean  of  the  colatitudes  of  the  two  points. 

If  only  one  crossing  of  a  boundary  is  found,  or  if  no  crossings  are  found  (or  SSJ/4  data  is 
missing),  the  boundaries  are  determined  as  follows.  Equatorward  boundaries  are  based  on  an 
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analytic  representation  of  the  Gusserihoven  et  al.  [1983]  boundary.  The  boundary  itself  is  a  circle 
whose  center  is  displaced  from  the  magnetic  pole.  The  radius  (in  degrees)  of  the  circle  is 


6](P1)  =  20.9°  +  1.70/> 

where  Px  is  a  "precipitation  index".  The  center  of  the  circle  is  located  at 


M^)  =  87.30-€.267°P] 

<p  cl(/»)  =  39-5°- 1.25°  Px  +  0.076° /f 

where  Xcl  is  magnetic  latitude  and<pcl  is  magnetic  local  time  (in  degrees). 

Poleward  boundaries  have  almost  the  same  form,  except  that  they  are  parameterized  in 
terms  of  a  separate  "precipitation  index"  P2. 

02(/>2)  =  13.4°+  1.7°  P2 


j89.2°+0.267°/>2, 

|90.8°-0.267°P2, 


P2<  3 
P2>  3 


^2(^2) 


cpcl  +  180°,  P2<3 

.9  ci’  ^2  —  3 


If  there  is  no  data,  or  if  no  boundaries  are  detected  in  the  data,  then?!  -P2-Kp.  However,  if  a 

single  crossing  is  detected,  then  the  precipitation  index  is  chosen  so  that  the  above  boundary 
matches  the  crossing. 
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L0MNE  '  Real  v  ?;|||||  y;  ±XXX.XX  ;  ;'||1|.:  degrees  longitude  of  northeast  comer  of  1  per  record 

X?I:':..  'M::'"':/-  .y  llfr?1!/.- :  ’  :l '  :;:r; '  ;f?;;  user  grid  (positive  east)  v-y  ;  {  : 


Record  #  Element  Name  Data  Type  ASCII  Format  Units  Description  Repetition  Rate  Notes 

16  SPACNG  Character  “N”or“G'’  N/A  “N”  indicates  user  will  specify  1  per  file  .  6 

^’.number  of  grid  points;  f  ■  •.  ||  f-  ;|l!:Illk  •• 

“G”  indicates  user  will  specify  size 
_ ;• _ _ _ of  crid  intervals 
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If  present,  record  1 1  must  contain  a  valid  file  name  of  no  more  than  32  characters.  The  file  will  be  written  to  the  default  directory  (the  one 
containing  PRISM.EXE). 

This  record  is  present  only  if  USRGRD  (record  12)  is  present  and  contains  “Y” 

This  record  is  present  only  if  SPACNG  (record  16)  is  present  and  contains  “N”. 

This  record  is  present  only  if  SPACNG  (record  16)  is  present  and  contains  “G”. 


05 

& 

I 

x 

o 

<N 

CO 

c 

G 


2 

o 


o 
c 

O 

I 

G 

g 
gS 
> 


c4 

t-H 

o 


<D 

X 


o 

> 

a 


c/2  C 

g  x 
jz  j  ^ 
S  c 
VO  O 

Uh  *-» 
O  c/5 

8  s 

£  s 

H  3 
S  3 
O 

J  <L> 

o  £ 

U  M 

o 


’E  8 

O 
o 

.8 


C/5 

.s 

*3 


H  o 

5  £ 

O  CTn 

|  i—i 

H 

00  g 

o 
<u 


o 

X 

H— > 

o 

X) 

K 

o 

<u 

g 

o 


& 

o 

.8 

X3 

+-» 

3 

,G 

<4H 

<U 

"O 

<L> 


4> 

X 


<D 

X 

d 

Oh 

G 

.2 

<D 

> 

g 

o 

o 

1 

I 


t-H 

a> 

o 

cx 

Oh 

<l> 

sc 


W) 

c 

*o 

g 

o 

Ph 

V5 

8 

o 

o 


g 

a> 


<4H 

<D 

O 

o 

o 

a> 

<D 

(75 


cq 

§ 

caT 

<4-1 

o 

G 


G 

O 

+-» 

G 

<U 

C/5 

4> 

Oh 


>*  15 
3  3 

o  o 
^  O 


3 

o 

o 

c/5 


£ 
4-» 
05 

G 

s 

T3  <D 


<D 

X 


G 

O 

C/5 

<D 

C/5 

G 

S 

CO 


r2 

05 

£ 

o 

G 

G 

G 


1 
II  s 

H  3 

G  G 


<N 

O 

o 

05 

G 

r-^ 

<3 

a-> 

05 

3 

+-> 

G 

O 

G 

O 

05 

o- 

ao 

CN 

3 

3 

s 

o 

O 

3 

40 

40 

'G 

T3 

2 

3 

"O 

CO 

G 

CO 

a 

a> 

o 

d 

G 

r-H 

G 

E 

Gh 

40 

o 

O 

C-) 

'O' 

G 

G 

O 

G 

3) 

c 

G 

4J> 

05 

3 

o 

*3 

±3 

<U 

g 

c/5 

05 

Uh 

Oh 

o 

05 

+-> 

o 

o 

•  ^ 

T3 

G 

II 

5-h 

o 

2 

CO 

CO 

o 

IX 


G 

3 

G 

o 

o 

C/5 

& 

js 

3 

3 

3 

o 

X 

C/5 

00 

<N 

TO 

Vh 

o 

u 

<L> 

m 

o 

2 

< 

2 

U 


uo 


vo  r-  oo  ov 


<n  m  Tj“ 


ON 

a 

o 

c 

T3 

C 

ca  ^ 
On  VO 


^  P  ’E 
o  o  5 


G 

O 


> 

r — V 

On 

(N 

i 

r- 

r4 


G* 

G*  On 
On  On 
ON  *-H 


G 

on 

on 


v-  V 

|l* 

E  o 
£  3 

cx  <io 

o  CO 


CO 


VO 


*o 
On 

_  On 

1  gj 
>  ^2 


<N  ^ 


O 
m  cD 


g  -H  <N  CD 

TO  G  G  C 
^  O  O  O 
M  *jd  *X3  X 
G  o  o  o 

t  t 
X  o  o  o 
U  u  o  o 


47 


m 


m 

06 

IP 


a 

c 

»2 

1  as 

S  Sn 

Pa  ^ 

co  *2 
»  JS 
f?  o 
fe  era 

S 

P, 

e 


© 

S=3 

c 

© 

U 


£d 


< 

z, 

1 

P 

b 

£ 

cd 

2 

r2 

E 


<u 

Pi 

cT 

04 


© 


9-n 

<D 

P 


i 

03 

O 


XX© 


^  <N 
© 
X 


•S1 

o 

C/3 

© 

Q 


03 

03 

© 

O 

o 

© 


O 

g 
© 

03  r£ 

CO 

cd 

>  03 

£ 
t> 


:§  T3  g 

~  £i 


s 

o 

ss 

03 

0) 

-*-> 

cd 

o 

•  fH 

X 

C 


G 

O 

x 

CO 


CO 


> 

X 

< 


TJ 
© 
g 
© 

P 

o  > 


a 


<D 

£f 


X 


© 

bD 

© 


co 

W 

O 

U 

< 

Q 


co 

pf 


0) 

X 

© 

CL, 


03 

.£ 

*3 

-4— » 

c 

o 

Q 

c3 


03 

+-* 

a 

© 

O 

S 

*L-h 

© 

© 

Q 

© 

X 


& 

s  j 

£ 

^4-> 

c 
o 
o 

i-^ 

cd 

X 


o 

.§ 

"a 
© 

03 

£  U 
Sd  E 
o  *+* 
©  <u 
Da  x: 

03  *  ‘ 


03 

~  © 
O  *-> 
^  Cd 

o  c 

i  ’■§ 

§  § 

o  w 

©  o 
bo  X 
G  © 
i  & 

<£  | 

a>  o 

■fi  S{ 

S-> 

=  S 

s « 


_  o 

L-h  o 


03 

*C3 

< 

!< 

!< 

D 

Z 

lz 

iz 

jz 

S3 

H 

< 

& 

8 


<N 


JJ 

hg 

JJ 

.a 

u2 

,2 

s 

|  IG 

i« 

la 

l«a 

I E 

|E 

Lh 

1 

1  S-! 

1  *“H 

1  Jh 

i  , 

i  ^ 

© 

<D 

© 

(U 

© 

<L> 

P 

1  P 

P 

P 

P 

i  p 

P 

Y-ccj 

■  r—i 

h 

h 

1- 

n 

*  t-H 

bJQ 

a 


co 

co 


c 

’3  .  . 

c  H 

8  CO 

o  O 


03  *-> 

£-3 

©  Q 

PX 

03  -4-* 


co 

CO 

co" 

w 

H-t 

co 

co 


CO 


bO 

.£ 

’c 

*3 

ts  *-» 

£  c 
o  0) 
o  - 

bis 

O  <o 

W  pi 

o  w 

•X  © 

©  o 
sz  S 

^  ^  03 

03  rtj 

©  ^  So 
<x  co  c3 
*1G  |— )  X> 
o  ^  C3 

<D 

CL  r£ 

03  *£3  T3 


bo  a 

j  c.l 

bJ) 

G  P 

G  fcn 

'£  ‘o 

3  fcj 

c  2 

■a  ^ 

l*G  <U 
i.S 

I  G  G 

3 

§  s 

O  V 

1  -4—*  (U 

1  4_>  1— , 

i  G  O 

J  o  o 

j  §  w 

0  «j 

.  ?-4 

b  ^ 

!  &  ^ 

1  &  G 

2  S* 

|  o  G 

|  O  -g 

o 

[  1 

1  0  3 

<u  U-t 

1  P  Gb 

.S3 

r°  "2 

O  G 

©  3 

x  .td 

1  X  \£ 

x:  ^ 

*-*  '*— * 
03  ^3 

03I  ~  cd 

«i  «  r 

03!  ^  £ 

PJ  Sf?  ? 

itt  c  ^ 

o  C  rt 

©  flj 

Ph  JZ!  iS 

m  t3  T5 


< 

2 


(N 

I  P! 

I  <N 

|  PS 

|  <N 

1  PS 

I  PS 

i  p 

2 

1  a> 

1  O 

1  <D 

i  2 

1  2 

I  aj 

1  <D 

0 

1  O 

1  0 

f  O 

I  0 

1  0 

l  0 

1  O 

G 

1  G 

I  G 

1  G 

1  G 

1  G 

1  G 

1  G 

(U 

1  <D 

I  © 

1  © 

1  O 

1  <D 

1  (D 

1  2 

© 

|  O 

I  <u 

1  &> 

1  © 

i  P 

j  QO 

l  2 

03 

03 

j 

03 

j 

03 

j 

1“ 

| 

03 

| 

03 

03 

| 

L-, 

l 

1  ?-4 

1 

1  ?— < 

1 

1  *-i 

1 

1  5-h 

1 

I  ?-< 

1 

I  ^ 

1 

1  h 

2 

1  2 

1  2 

1  2 

1  2 

1  2 

1  2 

1  2 

0 

1 

1  3 

1  <-> 

1  0 

1  <■> 

1  2 

l  2 

cd 

cd 

cd 

cd 

Cd 

cd 

cd 

3 

£ 

% 

% 

% 

S3 

% 

1  ^ 

1  ®3 

S3 

jC 

S3 

JC 

X! 

X 

X 

X 

u 

ju 

!U 

U 

|o 

lu 

1° 

|U 

lOi 


•P  O  JD 

W  O  3 
CL  r*  cd 

03  rS  ^3 


< 


is 

1 

1 

1 

I 

ffi 

1 

Ip 

< 

[P 

< 

a, 

]5 

< 

Is 

l< 

H 

< 

la. 

i< 

p 

I 

l< 

J&3 

p 

D 

CL 

1  w 

w 

°5 

H 

&o 

lb 

P 

JS 

J 

a 

P5 

D 

!2 

J 

3 

D 

VO 


--  & 
Sa  o 

O  — 1 
<U  a) 

8?^ 
03  -4-^ 


< 

Z 


03 

<D 

03 

C3 

X 

cd 

cd 

T3 

G 

.2 

O 

s 

<D 

O 

© 


i  pq 

iso  ’S 

c  5 


.£ 

'3 

G 

O 

o 

>> 

u 

o 

•w 

o 

© 

03 

D 

s 

p 

03 


< 


cd 

bX) 

s 

w 

j 

E 

>» 

X) 

-o 

T3 

© 

© 

G 


O 

"O 

o 

£ 

& 

ai 


© 


bJO 

G 

G 

’3 

^— < 

G 

O 

°  G 
>>  D 
Vh  — 


03 


2  ,3 

o 

(U 


o 
g 

^t-t 

■■§  s 

03  OC 

<U  OC 


o  cq 

© 

P  LS 

03  X 


00 


< 

2 


ON 


48 


e  cd 


eft 

eft  7P 
V  £ 
O  © 

o  ^ 

*  S 


£ 


O  eft 

ZO 
Q 

Is 

L  A  \ 


13 


Cft 


o  — <  jz 
>  o>  £ 
^  £  o 


<D 


<1> 


O 

S-H  £- 


cd 


O 


eft  £ 

•T— ( 


J8 


cs  cn  ^t 


1 

03 

.5 

S 

<r> 

a 

o 

£ 

*8 


o 

rf 

0> 


On  io 

^  ON 

5  2 

if 

6  i 

CO  h-> 
<N  ON 


CD 


49 


Version  1.5  30  May  1995  (No  changes  except  version  number) 
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names  are  not  always  descriptive  of  the  actual  contents  of  the  data  element  (for  historical  reasons). 

Because  PRISM  uses  FORTRAN  list  directed  input  for  all  numeric  input,  the  numeric  format  is  flexible.  The  formats  indicated  in  the  table  are 
suggested  formats  only.  Plus  (+)  signs  are  always  optional.  The  variables  may  be  separated  by  (1)  one  or  more  spaces,  (2)  a  comma  (,),  or  (3)  a 
comma  and  one  or  more  spaces. 


All  records  have  the  same  format.  PRISM  recognizes  the  end-of-file  (eof)  mark  as  an  “end-of-data”  indicator.  Thus  the  number  of  records, 
does  not  have  to  be  specified. 

PRISM  versions  through  1.4  ignore  Fi  layer  parameters.  However,  future  versions  way  make  use  of  these  parameters  so  space  is  reserved  in 
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All  records  have  the  same  format.  PRISM  recognizes  the  end-of-file  (eof)  mark  as  an  “end-of-data”  indicator. 

The  precise  definitbn  of  the  IIP  used  to  calculate  the  vertical  equivalent  TEC  is  not  critical  as  long  as  it  is  used  consistently. 
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Record  #  Element  Name  Data  Type  ASCII  Format  Units  Description  Repetition  Rate  Notes 

2  (“SSJ/4  YEAR  Integer  XXXX  years  4  digit  year  of  observation  1  per  record  2,3 

DATA”)  DAY  Integer  XXX  days  day  of  year  of  observation  1  per  record 

UT  Integer  XXXXX  seconds  Universal  Time  of  observation  1  per  record 

GLAT  Real  ±XX.XX  degrees  geographic  latitude  of  spacecraft  1  per  record  4 
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A  DMSP  input  file  for  PRISM  1.4  may  contain  three  kinds  of  data,  each  idenfified  by  a  unique  header  string: 
“SSIES  ION  DRIFT”  identifies  Drift  Meter  (DM)  data 
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Record  # _ Element  Name  Data  Type  ASCII  Format  Units  Description  Repetition  Rate  Notes 

9  LATO  Real  XXXX.XX  degrees  starting  latitude  of  grid  (southwest  1  per  file 

corner) 

LATF  Real  XXXX.XX  degrees  ending  latitude  of  grid  (northeast  1  per  file 

comer) 
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The  structure  of  the  PRISM  gridded  output  file  after  Record  10  depends  on  the  content  of  Record  10. 
For  IOUT=0,  Records  1 1-13  are  repeated  until  the  data  for  every  grid  point  has  been  written. 

For  IOUT=l  or  2,  Records  1 1  and  12  are  written  only  once. 

This  Record  is  repeated  until  the  contents  of  the  Data  Element  have  been  exhausted. 
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The  structure  of  the  PRISM  station  output  file  after  Record  2  depends  on  the  content  of  Record  2. 
For  IOUTS=0,  Record  4  is  repeated  NRECL  times. 
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